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--^ , The virialization of matter overdensities within dark energy sub- 

q-( systems is considered under a number of restrictive assumptions, namely 

(i) spherical-symmetric density profiles, (ii) time-independent quintessence 
bi)' equation of state parameter, w, and (hi) nothing but gravitational in- 

teraction between dark energy scalar field and matter. In addition, 
the quintessence subsystem is conceived as made of "particles" whose 
mutual interaction has intensity equal to G(l + 3w) and scales as the 
inverse square of their distance. The related expression of the self and 
tidal potential energy and formulation of the virial theorem for sub- 
systems, are found to be consistent with their matter counterparts, 
passing from —1 < w < —1/3 to w = 0. In the special case of fully 
clustered quintessence, energy conservation is assumed with regard to 
either the whole system (global energy conservation), or to the matter 
subsystem within the tidal potential induced by the quintessence sub- 
system (partial energy conservation). Further investigation is devoted 
to a few special cases, namely a limiting situation, w = —1/3, and 
three lower values, w = —1/2,-2/3,-1, where the last one mimics 
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the effect of a cosmological constant. The special case of fully clus- 
tered (i.e. collapsing together with the matter) quintessence is studied 
in detail, using a similar procedure as in Maor & Lahav (2005). The 
general case of partially clustered quintessence is considered in terms 
of a degree of quintessence de-clustering, £, < ( < 1, ranging from 
fully clustered (£ = 0) to completely de-clustered (£ = 1) quintessence, 
respectively. The special case of unclustered (i.e. remaining homoge- 
neous) quintessence is also discussed. The trend exhibited by the frac- 
tional (virialization to turnaround) radius, rj, as a function of the (i) 
fractional (quintessence to matter) "mass" ratio at turnaround, m, (ii) 
degree of quintessence de-clustering, (, and (iii) quintessence equation 
of state parameter, w, is found to be different from its counterparts 
reported in earlier attempts. In particular, no clear dichotomy with 
respect to the limiting situation of vanishing quintessence, r\ = 1/2, 
is shown when global or partial energy conservation hold in the spe- 
cial case of fully clustered quintessence, with rj > 1/2 preferred. The 
reasons of the above mentioned discrepancy are recognized as owing 
to (i) different formulations of the virial theorem for subsystems, and 
(ii) different descriptions of de-clustered quintessence, with respect to 
the reference case of fully clustered quintessence. 
keywords - Dark matter; Dark energy; 95.35.+d; 95.36.+X. 

1 Introduction 

Recent observations from anisotropies in the cosmic microwave background, 
large-scale structure surveys, Hubble parameter determinations, and Type 
la supernova results, allow narrow ranges for the values of cosmological 
quantities (e.g., Sievers et al. 2003; Rubigho-Martin et al. 2003; Spergel 
et al. 2003). The related "corcondance" cosmological model is consistent 
with a bottom-up picture (hierarchical clustering) of dark matter haloes, 
where smaller systems formed first from initial density perturbations and 
then merged with each other to become larger systems, or were tidally dis- 
rupted and accreted from successively formed, more massive neighbours. A 
main feature of the concordance model is that the dominant (about 70%) 
contribution to the present-day energy budget is a component with equation 
of state, p = wc 2 p, called dark energy, where p is the pressure, p the energy 
density, c the speed of the light, and w a dimensionless parameter which is, 
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in general, time-dependent. 

An obvious candidate, a scalar field, has to be light enough to vary slowly 
during a Hubble time, in such a way its potential energy can drive an ac- 
celerated expansion, just like during inflation. The varying field equation of 
state can then be tuned to lie in the observed range, and the related scalar 
field is sometimes called "quintessence", which may be conceived as a fifth 
component of the cosmic fluid after ordinary matter (baryons, leptons, and 
radiation) and nonbaryonic dark matter. The scalar field density fraction, 
Q q , can be made to decrease rapidly in the past, so as to pass easily the lens- 
ing constraints and to avoid discrepancies in the primordial nucleosynthesis 
abundances. For further details, see e.g., Amendola (2000). 

The parameter of the equation of state, w, depends on how the scalar 
field is slowly rolling down its potential. In the limit of a completely flat 
potential, w — — 1, the quintessence behaves as a cosmo logical constant, A > 
(Wetterich 1988; Peebles & Ratra 1988; Ratra & Peebles 1988). A subclass 
of quintessence models with constant w was proposed by Caldwell et al. 
(1998a,b) within the range, — 1 < w < —1/3. In fact, w > —1/3 corresponds 
to decelerate expansion while w < — 1 implies a phisical interpretation which 
is not still clear up today (e.g., Maor & Lahav 2005, hereafter quoted as 
ML05). Data on large-scale structures suggest a preferential range, — 1 < 
w < —0.6 (e.g., Weinberg & Kamionkowski 2003). Many other forms have 
been proposed for the shape of the potential of the scalar field, leading to an 
equation of state parameter that is dependent on the scale factor (see Peebles 
& Ratra 2003, for a review). For further details see e.g., Percival (2005). 

Dark energy affects not only the expansion rate of the background and the 
distance-redshift relation, but also the growth of structure. The formation 
rate of haloes, their evolution and their final characteristics are modified. 
Dark energy is therefore expected to have an impact on observables such as 
cluster number counts and lensing statistics due to intervening concentrations 
of mass on the line of sight of background sources. For further details see e.g., 
Horellou & Berge (2005). The assumption of an Einstein-de Sitter cosmology, 
Q m + Qg = 1, according to the concordance model, allows a much simpler 
analysis, but additional assumptions are needed. 

First, overdensities are conceived as spherical-symmetric, as initially done 
by Gunn & Gott (1972). The related top hat spherical collapse formalism 
has been proven to be a powerful tool for understanding the formation and 
the evolution of bound systems in the universe. For reasons of simplicity, 
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the special case of homogeneous overdensities is currently considered (e.g., 
Lahav et al. 1991; Wang & Steinhardt 1998; Iliev & Shapiro 2001; Horellou 
& Berge 2005; Percival 2005; ML05). 

Second, the quintessence equation of state is restricted to constant val- 
ues of the parameter, w. A time-dependent parameter, w, can result from 
a changing ratio of quintessence kinetic to potential energy which, in turn, 
is owing to the evolution of the scalar field potential (e.g., Caldwell et al. 
1998b). Though time- varying equations of state have been used in the lit- 
erature (e.g., Wetterich 1995; Amendola 2000; Battye & Weller 2003; Mota 
& van de Bruck 2004; Percival 2005), the special case of constant w makes 
considerable simplification (e.g., Caldwell et al. 1998b; Wang & Steinhardt 
1998; Weinberg & Kamionkowski 2003; Horellou & Berge 2005; ML05). 

Third, any couplings of quintessence to other fields are supposed to be 
negligibly small, so that the scalar field interacts with other matter only 
gravitationally. Different interactions can be taken into consideration (e.g., 
Wetterich 1995; Amendola 2000) at the price of less simple analysis involving 
a large number of parameters. 

In this view, the quintessence is conceived as lying between two limiting 
cases, namely (i) full clustering i.e. the scalar field responds to the infall in the 
same way as matter, and (ii) unclustering i.e. its sole effect is a tidal potential 
acting on the matter overdensity. The general case of partial clustering can 
be taken into consideration, at the expense of a much more complicated con- 
tinuity equation for the quintessence overdensity (e.g., Mota & van de Bruck 
2004; ML05). Even though it has been shown that quintessence cannot be 
perfectly smooth (e.g., Caldwell et al. 1998a,b), clustering is usually assumed 
to be negligible on scales less than about 100 Mpc (e.g., Wang & Steinhardt 
1998; Weinberg & Kamionkowski 2003; Battye & Weller 2003; Horellou & 
Berge 2005). It is therefore common practice to keep the quintessence ho- 
mogeneous during the evolution of the system. The effects of relaxing this 
assumption were explored in recent attempts (e.g., Mota & van de Bruck 
2004; Percival 2005). The real situation may safely be expected to lie be- 
tween two limiting cases, where the quintessence is clustering together with 
the matter and remains homogeneous, respectively. 

In spite of the simplifying assumptions mentioned above, still different 
authors use different expressions for (i) the gravitational potential induced 
by the quintessence, (ii) the self and the tidal potential energy, and (iii) the 
formulation of the virial theorem, which yields different results. On the other 
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hand, considerable effort has been devoted to the fact, that the quintessence 
is indistinguishable from a cosmological constant in the special case, w — — 1 
(e.g., Wang & Steinhardt 1998; Mota & van de Bruck 2004; Horellou & Berge 
2005; Percival 2005; but see also ML05; Wang 2006), but little consideration 
has been taken on the fact, that the quintessence is indistinguishable from 
matter in the special case, w = (e.g., Caldwell et al. 1988b; Iliev & Shapiro 
2001). Accordingly, the expression of the gravitational potential induced 
by the quintessence subsystem within a density perturbation, is expected to 
reduce to its matter counterpart as w — > 0. In fact, quintessence and matter 
gravitational potential are sometimes unified in a single expression, where 
w < —1/3 corresponds to the former case and w = to the latter (e.g., 
ML05). 

The current paper is restricted to the virialization of matter overdensi- 
ties within quintessence subsystems, under the simplifying assumptions men- 
tioned above and, in addition, the boundary condition that the quintessence 
potential reduces to the gravitational potential, when the pressure term is 
suppressed. The gravitational potential and potential energy terms for mat- 
ter and quintessence subsystem, are expressed in Sect. 2. The virial theorem 
for 2-component systems is specified in Sect. 3. The special case of fully clus- 
tered quintessence, and its extension to the general case of partially clustered 
quintessence, are investigated in Sects. 4 and 5, respectively. The related dis- 
cussion is performed in Sect. 6, and some concluding remarks are reported in 
Sect. 7. Further investigation on a few special arguments is developed in the 
Appendix. 

2 Matter and quintessence potential-energy 
terms 

Let us take into consideration homogeneous and spherical-symmetric over- 
densities made of both matter and quintessence. Let the quintessence equa- 
tion of state: 

p q = wc 2 p q ; (1) 

be restricted to constant values of the parameter, w, being p the pressure, 
p the density, c the light velocity in vacuum, and the index, q, denoting 
quintessence. Let any couplings of quintessence to other fields be negligibly 
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small, so that the scalar field interacts with other matter only gravitationally. 

The overdensity may be conceived as a two-component system made of 
matter and quintessence, respectively. The related gravitational potential is 
(e.g., Mota & van de Bruck 2004; ML05) 1 : 



where w m = 0, w q = w, R is the radius of the overdense sphere, and the 
index, to, denotes matter. 

With regard to matter component, Eq. (2) reduces to the well known 
expression of the gravitational potential induced by a homogeneous sphere 
on interior points, with the boundary condition that the maximum value is 
attained at the centre (e.g., MacMillan 1930, Chap. II, §29). In fact, given 
an interior point placed at a distance, r, from the centre, the potential in- 
duced from the sphere of radius, r, is y( int )(r) = (4n/3)Gpr 2 , according to 
MacLaurin's theorem, and the potential induced from the corona of radii, r 
and R, is V^ cxt \r) = 2irGp(R 2 — r 2 ), according to Newton's theorem, and 
the potential induced by the sphere of radius, R, at the point under con- 
sideration, is V(r) = y( mt )(r) + y( cxt )(r), according to Eq. (2) particularized 
to the matter subsystem. For further details see e.g., Caimmi (2003). An 
alternative expression is related to the boundary condition of a null potential 
at the centre (e.g., Iliev & Shapiro 2001; Horellou & Berge 2005). 

In general, the gravitational potential may be conceived as induced from a 
matter distribution where any two particles, idealized as mass points, interact 
with strenght, G, according to Newton's law, F m . m . = Grriimj/r 2 ^ and the 
dependence of the resulting force on the distance, dV/dr oc r, deduced from 
Eq. (2), is owing to the selected density profile, p = const. 

In this view, the gravitational potential induced by the quintessence may 
also be conceived as arising from a distribution where any two "particles", 
idealized as "mass points", interact with strenght, (1 + 3w)G, according 
to a Newton-like law, F mq . mq] = (1 + 3w)Gm qi m q j / rf^ where m q is the 
"quintessence mass". Then the results related to two-component matter 
distributions (e.g., Limber 1959; Brosche et al. 1983; Caimmi et al. 1984; 
Caimmi & Secco 1992) may be generalized to the case, where a subsystem is 
made of quintessence. 

1 A different sign convention is adopted here. 
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The self potential energy reads (e.g., Chandrasekhar 1969, Chap. 2, § 10): 

^u = -\f Pu(x 1 ,x 2 ,x 3 )V u (x 1 ,x 2 ,x 3 )d 3 S ; u = m,q ; (3) 
^ J s u 

where S u is the volume of u subsystem. 

The interaction potential energy, the tidal potential energy, and the resid- 
ual potential energy read (e.g., Caimmi & Secco 1992): 



W uv = ~( p u (x l ,x 2 ,x 3 )V v {x l ,x 2 ,x 3 )d 3 S ; (4) 
^ J s u 

f 3 dV 

V uv = p u (x 1 ,x 2 ,x 3 )J2 x rTr- !L( i 3 S ; (5) 



2JsJ 

p u {X 1 ,X 2 ,X 3 ) }_^X r 

Quv = V uv - W uv ; u = m,q ; v = q,m ; (6) 

respectively. 

The combination of Eqs. (2) and (4) yields: 

W uv W vu 

— — — = — — — ; u = m,q ; v = q,m ; (7) 

1 + 6W V 1 + 6W U 

in the special case of subsystems obeying the same equation of state, w u — w v , 
the interaction potential energy is symmetric with respect to the exchange of 
indices (e.g., MacMillan 1930, Chap. Ill, §76; Caimmi & Secco 1992; Caimmi 
2003). 

In the special case of concentric, spherical-symmetric, homogeneous sub- 
systems, the gravitational potential is expressed by Eq. (2), and Eqs. (3), (4), 
(5), and (7) reduce to (e.g., Caimmi 2003): 

16 3 CM 2 

n u = --7r 2 (l + 3w u )GplRl = --(l + 3w u )^ ; u = m, q ; (8) 

3 . . GMf m /5 2 1\ 

^ = --(1 + 3^)-^——^-?/ - -J ; u = m,q ; v = q,m ; (9) 



^, = --(1 + 3^,)^- ; (10) 



where the "mass" of the quintessence subsystem, M q , and the dimensionless 
ratios, m and y, are defined as: 

4<7T 

M q = —R qPq ; (11) 



and the indices, % and j, denote the inner and the outer component, respec- 
tively. 

Using the virial theorem for both subsystems and the whole system, to- 
gether with Eqs. (7), (9) and (10), the tidal potential energy, Vj i: takes the 
explicit expression: 



Vji = 



3 GMf m 
5 



Ri 



V 



x 



3 

T 2 
5y 2 -5 



1+ i w v-3 



5y 2 -l 



+ 2^-3 



(13) 



and a formal demonstration is provided in the next section. The combination 
of Eqs. (10) and (13) yields: 



V ji = V ij 



1 + 3wi 
1 + 3vjj 



(14) 



which makes a relation between tidal potential energies. 

Under the further restriction that the two subsystems are bounded by a 
single sphere, Ri = Rj = R or y = 1, Eqs. (9), (10) and (14) reduce to: 



mq 



w, 



qm 



1 + 3w g 1 + 3w r , 



V, 



mq 



V, 



qm 



1 + 3w a 1 + 3w r , 



3 GM m M q 
5 R 
3 GM m Mg 
5 R 



(15) 
(16) 



and the interaction potential energy coincides with the tidal potential energy, 
which is the limiting case currently used in the literature (e.g., Lahav et al. 
1991; Wang & Steinhardt 1998; Eliev & Shapiro 2001; Horellou & Berge 
2005; Percival 2005; ML05). 

Also it is worth of note (Percival 2005) that there is some confusion in the 
literature about the exact form of V mq = W mq , and the (1 + 3w q ) term has 
sometimes been neglected in the past, although it is included im more recent 
work (Battye & Weller 2004; Horellou & Berge 2005; Percival 2005; ML05). 
Owing to a different choice of the potential, the current expression of V mq 
coincides with its counterpart calculated in ML05, while the corresponding 
result appearing elsewhere (e.g., Battye & Weller 2004; Horellou & Berge 
2005; Percival 2005) is different by a factor, -1/2. 
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3 The virial theorem for matter and quintessence 
two-component systems 

It is worth recalling that the potential induced by the quintessence has been 
interpreted (Sect. 2) in terms of an interaction of strenght, (1 + w q )G, which 
depends on the inverse square distance. Accordingly, the virial theorem for 
the whole system reads (e.g., Landau & Lifchitz 1966, Chap. II, § 10): 

2T + Vl = ; (17) 

where f2 is the potential energy (e.g., MacMillan 1930, Chap. Ill, § 76; Caimmi 
& Secco 1992): 

tt = tt rn + W rnq + W qrn + tt q ; (18) 

and T is the kinetic energy. 

On the other hand, the virial theorem for subsystems reads (e.g., Limber 
1959; Brosche et al. 1983; Caimmi et al. 1984; Caimmi & Secco 1992; 
Caimmi 2003): 

2T U + Q U + V UV = ; u = m,q ; v = q,m ; (19) 

where T u is the kinetic energy of u-th subsystem. 

At this stage, let us assume that the quintessence subsystem can retain 
some form of kinetic energy, in such a way it is allowed to virialize by itself, 
within the tidal potential induced by the matter subsystem. In this view, 
the substitution of Eq. (18) into (17) yields: 

2T m + 2T q + Q m + Q q + W mq + W qm = ; (20) 

and the summation of the left-side member of Eq. (19) with its counterpart 
where the indices, u and v, are interchanged, produces: 

2^m ~\~ 2Tq -\- ~t~ ~t~ ^mq ~t~ V qm , (^-0 

finally, the combination of Eqs. (20) and (21) yields: 

V mq + V qm = W mq + W qm ; (22) 
or, without loss of generality: 

V uv = W uv + Q uv ; u = m,q ; v = q,m ; (23) 

Qmq Qqm i (^^) 
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where Q uv = V uv — W uv is the residual potential energy (e.g., Caimmi & 
Secco 1992). 

In the special case of concentric, spherical-symmetric, homogeneous sub- 
systems, the combination of Eqs. (9), (10), and (23) yields: 

3„ MMfm /5 5 2 \ 

^ = -5(1 + ^)^(4-^) ; (25) 

and the combination of Eqs. (23) and (24) produces: 

V ji = W ji + Q ji = W ji -Q ij ; (26) 

finally, using Eqs. (9) and (25) makes Eq. (26) coincide with (13). The further 
restriction that the two subsystems are bounded by a single sphere, Ri = 
Rj = R or y = 1, implies Q uv = and then V uv = W uv . 

The virial theorem, as expressed by Eqs. (17) and (19), is different from 
2T + Q m + 2V mq = currently used in the literature (e.g., Horellou & Berge 
2005; Percival 2005), due to the subtraction of Q q or the addition of V mq , 
respectively. On the other hand, the formulation 2T — Rd(Q m + V mq + V qm + 
Qq) I dR = used in ML05 is valid for the whole system, but cannot be used 
for subsystems. For further details, see Appendix A. 

In what follows, it shall be intended that the two fluids are made of 
matter (w m = 0) and quintessence (w q = w = const), and fill the same 
volume (R m — R q — R, y — 1). 



4 Fully clustered quintessence 

In the case of fully clustered quintessence, the quintessence field responds 
to the infall in the same way as matter, and the related continuity equation 
reads (ML05): 

p q + 3(l + w)-p q = ; (27) 

r 

where r is the radial coordinate. An integration from turnaround (r = -R max ) 
to a generic configuration r = R, using Eq. (11) yields: 



; (28) 
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which shows that the quintessence mass is decreasing with radius, and 
liniR_ > o M q (R) = 0, in the case under consideration, — 1 < w < —1/3. 

Owing to Eqs. (8), (15), and (18), the potential energy of the matter 
subsystem within the tidal field induced by the quintessence subststem, is: 



fim + W mq 
Mq 

mqm ~ M m 



|^[l + (l + 3 W )m gm ] 



(29a) 
(29b) 



where the mass ratio, m qm , changes with time, according to Eq. (28). At 
turnaround, R = R ma , x , the assumption of homogeneity makes total energy 
coincide with potential energy i.e. null kinetic energy (e.g., Horellou & Berge 
2005; Percival 2005; ML05). 

Owing to Eqs. (8), (15), (16), (19), and (28), the potential energy of the 
matter subsystem within the tidal field induced by the quintessence subsys- 
tem at virialization, is: 



(fir. 



7/ 



V, 



mq J vir 



3 GM 2 m 

5 Rrm ax 



1 (1 + 3w)m 
V 



l+3w 



Rvir 
Rmpt: 



m 



V 

M q (R 



Mr, 



(30a) 
(30b) 



where the indices, max and vir, denote turnaround and virialization, respec- 
tively. 

With regard to the quintessence subsystem, the counterpart of Eqs. (29a), 
and (30a), via (8), (15), (16), (18), (29b), and (30b) read: 



n + w t 



qm 



(n q + v qm \ 



5 R 
3GMH 



1 + 3w)m qm + m qm 
(1 + 3w)m 2 . m 



5 R n 



rj 



l+6w 



+ 



7] 



l+3w 



(31) 
(32) 



where the assumption of homogeneity makes total energy coincide with po- 
tential energy at turnaround. 

Using Eqs. (18), (29), and (31), the potential energy of the system at 
turnaround is: 



fimax 



-^-^[l + {2 + 3w)m+{l + 3w)m 2 ] ; 

O ilmfly 



(33) 
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according to ML05. 

Using Eqs. (15), (16), (18), (30), and (32), the potential energy of the 
system at virialization is: 



J_GMl 

i 'vir 



5 R n 



(2 + 3w)m (1 + 3w)m 2 



7] 



l+3w 



+ 



7] 



l+6w 



which is equivalent to: 



3GM4 1 

[l + mr]~ 



-3w\ 



5 IL 



rj 



) 1 + (1 + 3w)mrj' 



-3w 



(34) 



(35) 



in terms of radius and quintessence mass at turnaround. 

Keeping in mind that, in the case under discussion, the total energy equals 
the potential energy at turnaround and one half the potential energy at viri- 
alization, the requirement of energy conservation (e.g., Wang & Steinhardt 
1998; Weinberg & Kamionkonski 2003; Battye & Weller 2003; Horellou & 
Berge 2005; ML05; but see also Percival 2005) after combination of Eqs. (33) 
and (35) yields: 



rj = 



1 1+mr] 3w l + (l + 3w)mr] ; 
~2 1 + m l + (l + 3w)m 



(36) 



where the fractional radius, rj = -R v ir/-Rmax, appears on both the left-hand 
and the right-hand side. The solution of the transcendental Eq. (36) allows 
the knowledge of the virialized configuration. 

The counterpart of Eq. (36) in ML05, Eq. (16) therein, is different for 
the following reason. The expression, RdQ v [ T /dR, has been calculated at 
constant quintessence mass for an assigned configuration in the current ap- 
proach, while the quintessence mass, M q , has been derived with respect to 
the radius, R, in ML05. 

The fractional radius, rj = R w - U / R max , is necessarily restricted within the 
range, < rj < 1. In fact, distances cannot be negative and, on the other 
hand, r/ > 1 would contradict the definition of turnaround radius, R max . To 
gain more insight, let us write Eq. (36) under the form: 



_ 1 1 1 + x 1 + (l + 3w)x 

V - 2-/W - 2 i + m 1 + (l + 3w)m ' 

x = mi]~ 3w ; < x < m ; < f(x) < 2 ; 



(37a) 
(37b) 
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where the function, f(x), is studied in Appendix B. In the limiting case of 
a vanishing quintessence field, m — > 0, x — > 0, and Eq. (37a) reduces to 
7] = 1/2, which is the known result in absence of dark energy. 
The combination of Eqs. (37a) and (37b) yields: 

ir'/W 1 l + x l + {l + 3w)x 

(38) 



,mj 2 1 + m 1 + (1 + 3w)m ' 

and the virialized configuration is determined by the intersection point be- 
tween the curve on the left-hand and right-hand side, respectively. 

To gain further insight, let us take into consideration a few special cases, 
namely w — — 1, —2/3, —1/2, —1/3, the last to be thought of as an interesting 
limiting situation. The general case, — 1 < w < —1/3, is expected to show 
similar properties as in the closest among the above mentioned particular 
situations. 

In the special case, w = —1/3, Eq. (38) reduces to: 

x 1 1 + x 



m 2 1 + m 



(39) 



the solution of which is: 



= - = — U ; (40) 

m m + 2 

via Eq. (37b). The fractional radius, rj, for different values of the fractional 
mass, m, is represented in Fig. 1, top left. 

In the special case, w = —1/2, Eq. (38) reduces to: 

i\ 2 / 3 _ 1 2 + x-x 2 

m) ~ 22 + m-m 2 ' ^ ' 

the solution of which is determined by the intersection between curves on the 
left-hand and right-hand side. The fractional radius, 77, for different values of 
the fractional mass, m, is represented in Fig. 1, top right, where intersections 
are marked by open squares. 

In the special case, w = —2/3, Eq. (38) reduces to: 

x \ 1 / 2 1 1 — x 2 



rn 



) = 21^ • < 42) 



the solution of which is determined by the intersection between curves on 
the left-hand and right-hand side. The fractional radius, 77, for different 
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Figure 1: The fractional radius, rj = R v i r /R mSLX , for different values of the 
fractional mass, m = M g (i? max )/M m , with regard to global energy conser- 
vation and special choices of the equation of state parameter, w, from top 
left in the clockwise sense: —1/3, —1/2, —2/3, —1. The value of m related 
to each case can be read as the abscissa of the ending point of full curves, 
which occurs at F(x) = 1/2. The intersection of the two curves on the left- 
hand and right-hand side of Eq. (38), related to a virialized configuration, is 
marked by an open square. Curves corresponding to F(x) = [x /m)" 1 ^^ 
start from the origin. Curves corresponding to F(x) = (1/2) f(x) end at 
F(x) = 1/2, according to Eq. (38). The horisontal axis corresponds to dom- 
inant quintessence mass, m — > +oo. 
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values of the fractional mass, m, is represented in Fig. 1, bottom right, where 
intersections are marked by open squares. 

In the special case, w = — 1, Eq. (38) reduces to: 



the solution of which is determined by the intersection between curves on 
the left-hand and right-hand side. The fractional radius, 77, for different 
values of the fractional mass, m, is represented in Fig. 1, bottom left, where 
intersections are marked by open squares. 

An inspection of Fig. 1 shows the following features. 

(i) For assigned turnaround radius, the virialization radius is a decreasing 
function of the quintessence mass for w = —1/3, and an increasing 
function for w < —2/3. The trend is non monotonic for w = —1/2, 
with the occurrence of a minimum point. 

(ii) For assigned equation of state parameter, w, a virialized configuration is 
allowed for turnaround fractional mass within the range, < m < m , 
where m is the value for which the function on the right-hand side 
of Eq. (38) is divergent. In the limit, m — > m , density perturbations 
turn around and virialize at infinite radius, while density perturbations 
where m > mo cannot virialize. 

(iii) The threshold fractional mass ranges from m — > +00 (w = —1/3) to 
m = 2 (w = -1/2); m = 1 (w = -2/3); and m = 1/2 (w = -1). 
Accordingly, the minimum amount of quintessence within a density 
perturbation, necessary to prevent matter virialization, is an increasing 
function of the equation of state parameter, w. 

The above results rely on the assumption, that both energy conservation 
and virialization hold for the whole system (ML05). If, on the other hand, 
only the matter virializes i.e. the kinetic energy of the quintessence subsystem 
is null, the procedure should be repeated using Eqs. (29) and (30) instead of 
(33) and (34). The result is: 




(43) 




(44a) 



x = mi] 



—3w 



< x < m 



(44b) 
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where the virialized configuration is defined by the intersection of the curves 
on the left-hand and right-hand side of Eq. (44a), the latter being straight 
lines. 

The counterpart of Eq. (44a) in ML05, Eq. (17) therein, is different for the 
following reason. Though no indication is provided in ML05 on how Eq. (17) 
therein has been derived, it can be seen that it is sufficient to omit the terms 
U12 and U 2 2 (V mq and Q q in the current notation) in the expression of the 
potential energy, Eq. (4) therein. On the contrary, using Eqs. (29) and (30) 
implies the omission of V qm and Q q , which explains the different results. 

In the special case, w = —1/3, Eq. (44a) reduces to rj — 1/2, which is the 
result for matter universes. 

In the special case, w = —1/2, Eq. (44a) reduces to: 



the solution of which is determined by the intersection between curves on the 
left-hand and right-hand side. The fractional radius, 77, for different values of 
the fractional mass, m, is represented in Fig. 2, top right, where intersections 
are marked by open squares. 

In the special case, w = —2/3, Eq. (44a) reduces to: 



the solution of which is determined by the intersection between curves on 
the left-hand and right-hand side. The fractional radius, 77, for different 
values of the fractional mass, m, is represented in Fig. 2, bottom right, where 
intersections are marked by open squares. 

In the special case, w — — 1, Eq. (44a) reduces to: 



the solution of which is determined by the intersection between curves on 
the left-hand and right-hand side. The fractional radius, 77, for different 
values of the fractional mass, m, is represented in Fig. 2, bottom left, where 
intersections are marked by open squares. 




(45) 




(46) 




(47) 
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Figure 2: The fractional radius, rj = -Rvi r /-R max , for different values of the 
fractional mass, m = M q (R max )/M m , with regard to partial energy con- 
servation and special choices of the equation of state parameter, w, from 
top left in the clockwise sense: —1/3, —1/2, —2/3, —1. The value of 
m related to each case can be read as the abscissa of the ending point 
of full curves, which occurs at F(x) = 1/2. The intersection of the two 
curves on the left-hand and right-hand side of Eq. (44a), related to a virial- 
ized configuration, is marked by an open square. Curves corresponding to 
F(x) = (x/m)^ 1 ^ 3 ^ start from the origin. Straight lines corresponding to 
F(x) = (l/2)[l + (l + 3w)x]/[l + (l + 3w)m] end at F(x) = 1/2, according to 
Eq. (44a). The horisontal axis corresponds to dominant quintessence mass, 
m — > +oo. 



17 



An inspection of Fig. 2 shows similar features as in Fig. 1 (see above), with 
the following difference. For an assigned turnaround radius, the virialization 
radius is an increasing function of the quintessence mass for w < —1/3, and 
independent of the quintessence mass (as in matter universes) for w = —1/3. 
Accordingly, a monotonic trend is exhibited. 

The results plotted in Fig. 1 and Fig. 2 imply the assumption of energy 
conservation with regard to the whole system and to the matter subsystem, 
respectively, the latter requirement being more restrictive than the former. 
For a more detailed discussion on energy conservation in density perturba- 
tions with matter and dark energy, see ML05 and Percival (2005). 

The comparison of the current results (Figs. 1-2) with their counterparts 
in ML05 (Fig. 2 therein, related to w = —4/5) shows that the trend of the 
fractional radius, i], as a function of the fractional mass, m, is increasing in 
both cases with regard to the whole system, while is increasing in the current 
paper and decreasing in ML05 with regard to the matter subsystem. A closer 
agreement could occur within the range, — 1/2 < w < —1/3. The reasons of 
the above mentioned discrepancy are due to different assumptions made in 
the current paper and in ML05, as explained above. 



5 Partially clustered quintessence 

In the case of partially clustered quintessence, the quintessence field responds 
to the infall to a lesser extent with respect to the matter. In the limiting 
situation of unclustered quintessence, the related continuity equation reads 
(ML05): 

p q + 3(l+w)-p q = ; (48) 

(X 

where a is the cosmological scale factor, which is related to the redshift, z, 
via a oc (1 + z)' 1 . In the special case, w = —1, Eqs. (27) and (48) coincide, 
yielding p q = const or, in other terms, a cosmological constant. 

An integration from turnaround (r = R max ) to a generic configuration 
r = R, using Eq. (11) yields: 

M q (R) = M q (R max ) (-—) ; (49a) 

\ JX m ax ' \ ^max / 

— = 1 + Zmax ; (49b) 

"max I -\~ Z 
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which shows that the quintessence mass is decreasing with radius, and 
liniR_ > o M q (R) = 0, in the case under consideration, — 1 < w < —1/3. 

In the general case of partially clustered quintessence, it may safely be 
expected that the virialized configuration of matter subsystem lies between 
the limiting situation of fully clustered and unclustered quintessence, respec- 
tively. It is worth remembering that energy conservation can no longer be 
assumed, as in presence of unclustered quintessence (e.g., Horellou & Berge 
2005; Percival 2005; ML05). Without loss of generality, let us suppose that 
a specified virialized configuration is attained along the following steps: (i) 
from turnaround to virialization related to fully clustered quintessence, and 
(ii) from virialization related to fully clustered quintessence to virialization 
related to partially clustered quintessence. 

To perform the transition between the above mentioned steps, the quintessence 
mass must instantaneously change from M g (_R^ ir ) to M 9 (i? vir ), where the pres- 
ence and the absence of the prime denotes the virialized configuration related 
to fully or partially clustered quintessence, respectively. 

The changes in self, tidal, and interaction potential energy, keeping in 
mind Eqs. (8), (15), and (16), are: 

AO 3GMI . 

vir vir 

+ |(l + 3u,) GM "^' fl "'» ; (51) 
5 R vir 



AV — AW - Z GMmM^) 

O XL,,: 



| 3 GM m M q {R vir ) 

5 -^vir 

AQ q = --(l + 3w) ^ + -{l + 3w) ^ ; (53) 

O ^vir "^vir 

from which the changes in total energy within the matter subsystem and the 
whole system, can be deduced. 
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To this aim, let us define the fractional masses and the fractional radii: 

M q (R vlI ) M q (R' vh ) - M q (R vlI ) 

M q (R max ) ' ^ M q (R nmx ) ' l ° 4j 

^ = |^ ; A V = ^-^ ; (55) 
and the combination of Eqs. (28), (54), and (55) yields: 

A M <?(Kir) = ( Ku_\ 3W = -3w . ^ 

max ) V Rmax J 

r ] + Ar ] = ^=r ]FC ; (57) 

in the special case of fully clustered quintessence, A/i = 0, the results of 
Sect. 4 continue to hold provided 77 is replaced by i] FC therein. 

The substitution of Eqs. (30) and (54)-(57) into (50)-(53) produces: 

AQ m — ; (58) 
o CM 2 

AV mq = AW mq = --(l + 3w)—^mA f i ; (59) 

O -^vir 

AV qm = AW qm = --—^mA^ ; (60) 

r£ vir 
o CM 2 

AQ q = --(l + 3w)— ^m 2 [2/iA/i + (A/i) 2 ] ; (61) 

-O-vir 

in terms of virialized configurations related to fully clustered quintessence 
and to the change, A/z. 

At this stage, the virialized configuration of the matter subsystem can be 
determined using a similar procedure with respect to the special case of fully 
clustered quintessence. The change in total energy, related to the transition 
from virialized configurations where the quintessence is fully clustered (R = 
R' vir ) to virialized configurations where the quintessence is partially clustered 
(R = Ryir), reads: 

^n(R' vir ) + An = ±n(R vir ) ; (62a) 
AQ = AQ m + AV mq + AV qm + AQ q ; (62b) 
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and the combination of Eqs. (34), (54)-(57), and (58)-(62) yields: 

1 



5 -Rmax 

3 GMj 

5 Rry\ ax 



[2 + 3w)m (1 + 3w)m 2 



'/FC 



1 (2 + 3w)m/i (1 + 3to)m 2 /i 2 



Afi = — 



3GM 2 



5 R n 



[2 + 3w)m + (1 + 3w)m 2 (2/i + A//)] 



A// 

VFC 



(63) 
(64) 
(65) 



in terms of radius and quintessence mass at turnaround. 

Keeping in mind that, in the special case of fully clustered quintessence, 
the total energy equalizes the potential energy at turnaround and one half 
the potential energy at virialization, the substitution of Eqs. (33), (64), and 
(65) into (62a) produces: 



[1 + (2 + 3w)m + (1 + 3w)m 2 } + [(2 + 3w)m + (1 + 3w)m 2 (2rj^ w - A/i)] 

= hl + {2 + 3w)mfi + {l + 3w)m 2 fi 2 ]- ; 
2 rj 



'Vfc 
(66) 



where, using Eq. (56) and performing some algebra, the term within brackets 
on the right-hand side of Eq. (66) may be cast into the form: 



1 + (2 + 3w)mfi + (1 + 3w)m 2 fx 2 = 

1 + (2 + 3w)m%$° + (1 + 3w)m 2 rfcg° - 

<p(w,m,r] FC , A/i) = 

m[(2 + 3w) + 2(1 + 3w)mrfcg° - (1 + 3w)mAfi]Afi 

and the combination of Eqs. (66) and (67) yields: 



w,m,r] FC ,Afi) ; (67a) 
(67b) 



V= 2 



1 [1 + (2 + 3w)mri^ w + (1 + 3w)m 2 r]^ w } - 



[1 + (2 + 3w)m + (1 + 3w)m 2 } + rj^cj) 



(68) 



which depends on the parameters, w, m, rjpc, and A/i. 

In the special case of fully clustered quintessence, A/i = i.e. \x = 
i]pQ W , and Eq. (67) reduces to (36) where r) has to be replaced by t/fc- The 
solution of the transcendental Eq. (68) allows the knowledge of the virialized 
configuration for fixed A/i. 
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The counterpart of Eq. (68) in ML05, Eq. (23) therein, is different for 
reasons discussed in Sect. 4 and, in addition, due to a different choice of 
the parameter related to partial clustering, 7 therein instead of A/x or A77, 
according to Eq. (56). 

The fractional radius, 77 = R v h/R max , is necessarily restricted within the 
range, < 77 < 1. In fact, distances cannot be negative and, on the other 
hand, 77 > 1 would contradict the definition of turnaround radius, i? max - In 
addition, the particularization of Eqs. (28) and (49) to the related virialized 
configurations, shows that M g (i? vir ) < M q (R' vir .) within the range of interest, 
— 1 < w < —1/3, which implies A/i > via Eq. (54). Accordingly, < A/i < 
t]pq w , and the sign of is opposite to the sign of AQ via Eq. (65) which, in 
turn, is equal to the sign of 77 — t/fc- 

To gain more insight, let us express the increment, A/i, in terms of the 
variable, r] FC , and a degree of quintessence de-clustering, (, as: 

^ = Cvfc w ; < c < 1 ; (69) 

the substitution of Eq. (69) into (67b) yields: 

(f>(w, m, 77 FC , A/i) = (2 + 3w)m<:vFC W + (1 + 3w)m 2 C(2 - C)?7f"c" I (70) 

where the effect of quintessence partial clustering is expressed by the param- 
eter, (, where the limit of fully clustered and fully de-clustered quintessence 
corresponds to ( = and ( = 1, respectively. The latter, of course, has no 
physical meaning, as quintessence cannot be devoided from the volume filled 
by the matter subsystem. On the other hand, it makes an upper limit to the 
quintessence de-clustering parameter, (. Further details on the function, 0, 
are given in Appendix C. 

The above results may be summarized as follows. Given a density pertur- 
bation with assigned value of quintessence equation of state parameter, w, 
— l<w< —1/3, quintessence to matter mass ratio at turnaround, m, virial- 
ized to turnaround size ratio in the special case of fully clustered quintessence, 
^fc, < i]fc < 1, and degree of quintessence de-clustering, (, < ( < 1, the 
related virialized to turnaround size ratio in the case of partially clustered 
quintessence, 77, is expressed by Eqs. (68), (69), and (70). 

Let us repeat that the special case of fully de-clustered quintessence, ( = 
1, makes a convenient upper limit but, on the other hand, it has little physical 
meaning. A true upper limit is related to the special case of unclustered 
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quintessence, ( = £ unc , where the combination of Eqs. (49), (54), (55), (56), 
and (69) yields: 

ri = (i-U c ) 1/3 v^(^^) 1+W ; (71) 

which, in addition to Eq. (68), makes a further constraint for determining the 
virialized configuration, with regard to selected turnaround and virialization 
epoch. 

To get further insight, let us express the fractional radius, rj, in terms of 
the parameter, x = mrjp^ w . Accordingly, Eq. (68) reads: 

_ 1 [1 + (2 + 3w)x + (1 + 3w)x 2 } - 0(x; w, Q 

77 2 [1 + (2 + 3w)m + (1 + 3w)m 2 ] + (x/my/^^w, () ' 

where the function, <f>(x; w, (), is expressed by Eq. (114), Appendix C. 

The dependence, 77 = i](m), is represented in Fig. 3 with regard to same 
cases as in Fig. 1, but different degrees of quintessence de-clustering, ( = 
(squares), 0.25 (triangles), 0.5 (asterisks), 0.75 (crosses), and 1 (diamonds). 
The starting point on the left, (0,0.5), marks the limiting situation of vanish- 
ing quintessence, and necessarily coincides in all cases. For sufficiently high 
values of quintessence equation of state parameter, w < —1/3, lower curves 
correspond to less clustered quintessence and vice versa. For sufficiently low 
values of quintessence equation of state parameter, —1 < w < —2/3, lower 
curves correspond to more clustered quintessence and vice versa. For in- 
termediate values of quintessence equation of state parameter, w ~ —1/2, 
different curves intersect one with the other, but a minimum point occurs 
in the whole range between fully clustered (( = 0) and fully de-clustered 
(C = 1) quintessence. In any case, the general trend remains unchanged with 
respect to the limiting situation of fully clustered quintessence. 

The above results rely on the assumption, that both energy conservation 
and virialization hold for the whole system (ML05). If, on the other hand, 
only the matter virializes i.e. the kinetic energy of the quintessence system is 
null, the procedure should be repeated with regard to the matter subsystem 
only. Accordingly, Eqs. (62), (63), (64), and (65) are turned into: 

^m(-Rvir) + AE m = E m (R vir ) ; (73) 
AE m = An m + AV mq ; (74) 
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Figure 3: The fractional radius, i] = R v i T /R max , with regard to global energy 
conservation, for different values of the fractional mass, m = M q (R max ) / M m , 
the quintessence de-clustering parameter, (, and the quintessence equation of 
state parameter, w. From top left in the clockwise sense: w = —1/3, —1/2, 
—2/3, —1. Caption of symbols: squares, ( = (fully clustered quintessence, 
as in Fig. 1); triangles, ( = 1/4; asterisks, ( = 1/2; crosses, ( = 3/4; dia- 
monds, C = 1 (fully de-clustered quintessence). The value of m related to 
each symbol, starting from the left, is m = i/4, i — 0,1, 2, and the special 
case, i = 0, corresponds to a vanishing quintessence. 
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2E m (R' vil ) — — 



3 GM, 



5 -Rmax ^FC 



1 + 



3w)m 



„3ui 
^FC 



(75) 



2E m (R vir ) = --— -^-[1 + (1 + 3«;H ; (76) 
A^ = -»l + 3 W )^; (77) 

where E m = (Q m + W mq )/2 = (Q m + V mq )/2 is the total energy of the virial- 
ized matter subsystem, and AE m is the energy change due to the transition 
from a virialized configuration where the quintessence is fully clustered, to 
its counterpart where the quintessence is partially clustered. Using energy 
conservation in the former alternative, the substitution of Eqs. (75), (76), 
and (77) into (73) yields after some algebra: 

! [! + (! + 3w)mr]^ w } ~ <f> . (7g) 
V 2 [1 + (1 + 3w)m] + ?7fc0 ' 
0( W ,m,C) = (l + 3 W )mC^ ; (79) 

where Eq. (56) has also been used. 

In terms of the parameter, x = m-q^c" ■, Eqs. (78) and (79), translate into: 

„ = - [! + (! + 3*0*]-* . (8Q) 

' 2 [1 + (1 + 3w)m\ + (x/my/( 3w )<j) ' V ; 

<Kx,w,t) = (l + 3w)(x ; (81) 

in the limit of fully clustered quintessence, ( = 0, = 0, Eq. (80) coincides 
with Eq. (44a). 

The dependence, 77 = r](m), is represented in Fig. 4 with regard to the 
same cases as in Fig. 2, but different degrees of quintessence de-clustering, 
( = (squares), 0.25 (triangles), 0.5 (asterisks, 0.75 (crosses), and 1 (dia- 
monds). The starting point on the left, (0, 0.5), marks the limiting situation 
of vanishing quintessence, and necessarily coincides in all cases. It is ap- 
parent that lower curves correspond to more clustered quintessence and vice 
versa. In the limit, w — > —1/3, virialized configurations coincide with their 
counterparts in absence of quintessence, independent of the degree of clus- 
tering. On the other hand, the dependence is enhanced as the quintessence 
equation of state parameter, w, attains lower values up to —1. In any case, 
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Figure 4: The fractional radius, r\ = -Rvi r /-R ma x, f° r different values of the 
fractional mass, m = M q (R max )/M m , the quintessence de-clustering param- 
eter, (, and the quintessence equation of state parameter, w, with regard to 
partial energy conservation. Captions as in Fig. 3. 
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the general trend remains unchanged with respect to the special situation of 
fully clustered quintessence. 

In dealing with fully clustered quintessence, an inspection of Figs. 3 and 4 
shows that partial energy conservation yields larger virialized configurations, 
with respect to their counterparts where global energy conservation holds. 
An opposite trend is found in ML05, for w = —4/5, in the case of fully 
clustered (Fig. 2 therein) and unclustered (Fig. 4 therein) quintessence. The 
reasons of the above mentioned discrepancy are due to different assumptions 
made in the current paper and in ML05, as explained in Sect. 4. 

The dependence, r\ = 77(C), is represented in Fig. 5 with regard to both 
global (squares) and partial (diamonds) energy conservation (in the special 
case of fully clustered quintessence), form = 1/4 and w = —1/3, —1/2, —2/3, 
— 1. In both cases, larger fractional radii are attained for (algebraically) lower 
quintessence equation of state parameters, w. In addition, with regard to 
fully clustered quintessence, partial energy conservation yields larger virial- 
ized configurations when compared to their counterparts where global energy 
conservation holds. Accordingly, partial energy conservation implies a larger 
amount of kinetic energy in matter subsystems, with respect to global energy 
conservation. An opposite result has been found in ML05 (Fig. 3 therein) for 
w = —4/5 and m = 1/5. The reasons of the above mentioned discrepancy 
are due to different assumptions made in the current paper and in ML05, as 
explained in Sect. 4. 

The dependence, 77 = rj(w), is represented in Fig. 6 with regard to both 
global (squares) and partial (diamonds) energy conservation (in the special 
case of fully clustered quintessence), for m — 1/4 and ( = (fully clustered 
quintessence), 1/4, 3/4, 1 (fully de-clustered quintessence). In both cases, 
the trend remains unchanged as the degree of de-clustering gets increased. 
Other features have already been discussed above. A different result has been 
found in ML05 (Fig. 5 therein), for w = —4/5 and m = 1/5. The reasons of 
this discrepancy are due to different assumptions made in the current paper 
and in ML05, as explained in Sect. 4. 

6 Discussion. 

In absence of a full knowledge on the dark energy, an investigation on the viri- 
alization of spherical overdensities appears to be physically meaningful, even 
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Figure 5: The fractional radius, r\ = i? vir /i? max , as a function of the 
quintessence de-clustering parameter, (, for a fractional mass, m = 
M q (R max ) I M m , and values of the quintessence equation of state parameter, 
w, from top left in the clockwise sense: —1/3, —1/2, —2/3, —1. Global and 
partial energy conservation (in the special case of fully clustered quintessence) 
correspond to squares and diamonds, respectively. 
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Figure 6: The fractional radius, 77 = R v - n /R maiX , as a function of the 
quintessence equation of state parameter, w, for a fractional mass, m — 1/4, 
and values of the quintessence de-clustering parameter, (, from top left in 
the clockwise sense: (fully clustered quintessence), 1/4, 3/4, 1 (fully de- 
clustered quintessence). Global and partial energy conservation (in the spe- 
cial case of fully clustered quintessence) correspond to squares and diamonds, 
respectively. 
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if particular assumptions are made. More specifically, the potential induced 
by quintessence has been conceived as due to quintessence "particles" , inter- 
acting via an inverse square distance law, with strength equal to (f + 3w)G, 
— l<w< —1/3. Accordingly, the potential induced by a quintessence parti- 
cle scales as an inverse distance, V(r p ) oc l/r p , which makes the formulation 
of the virial theorem, 2T = kQ, V(r p ) oc r*, k = —1 (e.g., Landau & Lifchitz 
1966, Chap. II, §10). 

With regard to the transition from turnaround to virialization, the refer- 
ence case of energy conservation may be related to either the whole (quintessence 
+ matter) system, or the matter subsystem within the tidal field induced by 
the quintessence subsystem. According to ML05, global energy conservation 
occurs in the special case of fully clustered quintessence. The same is as- 
sumed with regard to partial energy conservation, even if it is not clearly 
specified in ML05 [Eq. (17) therein]. 

Bearing in mind the reference cases of global and partial energy conserv- 
tion, different degrees of quintessence clustering may be expressed via the 
de-clustering parameter, (, < ( < 1, which is related to the change in mass 
with respect to a selected reference case, according to Eqs. (54)-(57) and (69). 
The limiting cases, ( — 0, 1, represent fully clustered and fully de-clustered 
quintessence, respectively. 

Although typical values of the fractional mass at turnaround, m = 
M q (T max )/M m , do not exceed a value of about 0.9 (e.g., ML05), the analy- 
sis of a larger range shows interesting features, namely (i) the existence of a 
threshold in m, increasing with the quintessence equation of state parameter, 
w, above which no virialized configuration is allowed, see Figs. 3 and 4, and 
(ii) the occurrence of a non monotonic trend in the fractional radius, 77, as 
a function of the fractional mass, m, with the occurrence of a minimum, for 
w = —1/2, a decreasing monotonic trend for w = —1/3, and an increasing 
monotonic trend for w = —2/3,-1, with regard to global energy conservation 
in the limit of fully clustered quintessence, see Fig. 3. 

The dependence, 77 = 77(C), related torn — 1/4, is also monotonic for both 
global and partial energy conservation in the limiting case of fully clustered 
quintessence, see Fig. 5. Owing to Eq. (71), the special case of unclustered 
quintessence occurs within the whole range, < £ < 1, according if 1 < 
a vir /a max < 1500 and/or — 1 < w < —1/3. 

The dependence, 77 = f](w), shows a monotonic trend for both global and 
partial energy conservation in the limiting case of fully clustered quintessence, 
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see Fig. 6, with increasingly different trends as ( increases. It can be seen 
that virialized configurations cannot take place for both sufficiently low w 
and sufficiently large (. 

In any case, more extended virialized configurations occur for partial 
(with respect to global) energy conservation in the limiting case of fully 
clustered quintessence, see Figs. 5 and 6. Accordingly, the latter alternative 
yields bounder configurations with respect to the former. In addition, an in- 
creasing degree of quintessence de-clustering makes more extended virialized 
configurations provided partial energy conservation holds in the special case 
of fully clustered quintessence, see Fig. 5. When global energy conservation 
holds in the special case of fully clustered quintessence, a similar trend occurs 
for sufficiently low quintessence equation of state parameters, w < — 2/3. No 
change appears for w = —2/3, while the trend is reversed for w > —2/3. 

The above results are different from their analogon in ML05, even if the 
expression of the potential energy is the same as in the current paper. The 
discrepancy could be due to (i) different formulations of the virial theorem, 
as discussed in Sect. 4 and in Appendix A, with regard to both the whole sys- 
tem and the matter subsystem, and (ii) different descriptions of de-clustered 
quintessence, by defining a quintessence de-clustering parameter in the cur- 
rent paper, and modifying the quintessence continuity equation in ML05, 
Eqs. (9)-(ll) therein. 

If the fractional radius, i], can be deduced from observations, it (to- 
gether with the related uncertainty) translates into the plane of Figs. 3-6 
as a horizontal band, and the value of the fractional mass, m, the degree of 
quintessence de-clustering, (, and the quintessence equation of state param- 
eter, w, can be constrained. In addition, global and partial energy conserva- 
tion in the special case of fully clustered quintessence, can be tested. 

7 Conclusion 

The virialization of matter overdensities within dark energy subsystems has 
been considered under a number of simplifying assumptions, namely (i) 
spherical-symmetric density profiles, (ii) time-independent quintessence equa- 
tion of state parameter, w, and (iii) sole gravitational interaction between 
dark energy scalar field and matter. The quintessence subsystem has been 
conceived as made of "particles" whose mutual interaction has intensity equal 
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to (1 + 3w)G and scales as the inverse square of the distance. The related 
expression of self and tidal potential energy and formulation of the virial 
theorem for subsystems, have been found to be consistent with their matter 
counterparts, passing from — 1 < w < — 1/3 to w — 0. In the special case 
of fully clustered quintessence, following ML05, it has been assumed either 
global energy conservation related to the whole system, or partial energy 
conservation related to the matter subsystem within the tidal field induced 
by the quintessence subsystem. Further investigation has been devoted to a 
few special cases, namely a limiting situation, w = —1/3, and lower values, 
w = —1/2,-2/3,-1, where the last one mimics the presence of a cosmo log- 
ical constant. 

The special case of fully clustered (i.e. collapsing together with the mat- 
ter) quintessence has been studied in detail, following a similar procedure 
as in ML05. The general case of partially clustered quintessence has been 
considered in terms of a degree of quintessence de-clustering, (, < ( < 1, 
ranging from fully clustered (( = 0) to completely de-clustered (( = 1) 
quintessence, respectively. The special case of unclustered (i.e. remaining 
homogeneous) quintessence has also been discussed. 

The trend exhibited by the fractional radius, i], as a function of the 
fractional mass, m, the degree of quintessence de-clustering, (, and the 
quintessence equation of state parameter, w, has been found to be different 
from its counterpart reported in earlier attempts (e.g., ML05). In particu- 
lar, no clear dichotomy with respect to the limiting situation of a vanishing 
quintessence, r\ = 1/2, has been shown when global or partial energy conser- 
vation holds in the special case of fully clustered quintessence, with rj > 1/2 
preferred. The reasons of the above mentioned discrepancy have been rec- 
ognized as owing to (i) different formulations of the virial theorem, and (ii) 
different descriptions of de-clustered quintessence, with respect to the case 
of fully clustered quintessence. 

References 

[1] Amendola, L. 2000, PhRvD 62, 043511 

[2] Battye, R.A., Weller, J. 2003, PhRvD 68, 083506 

[3] Brosche, P., Caimmi, R., Secco, L. 1983, A&A 125, 338 



32 



[4] Caimmi, R. 2003, AN 324, 250 

[5] Caimmi, R., Secco, L., Brosche, P. 1984, A&A 139, 411 
[6] Caimmi, R., Secco, L. 1992, ApJ 395, 119 

[7] Caldwell, R.R., Dave, R., Steinhardt, P.J. 1988a, ApSS 261, 303 

[8] Caldwell, R.R., Dave, R., Steinhardt, P.J. 1988b, PhRvL 80, 1582 

[9] Chandrasekhar, S. 1969, Ellipsoidal Figures of Equilibrium, Yale Univer- 
sity Press, New Haven 

[10] Gunn, J.E., Gott, J.R., III 1972, ApJ 176, 1 

[11] Horellou, C, Berge, J. 2005, MNRAS 360, 1393 

[12] Iliev, I.T., Shapiro, PR. 2001, MNRAS 325, 468 

[13] Lahav, O., Lilje, P.B., Primack, J.R., Rees, M.J. 1991, MNRAS 251, 
128 

[14] Landau, L., Lifchitz, E. 1966 Mecanique, Mir, Moscow 
[15] Limber, D.N. 1959, ApJ 130, 414 

[16] MacMillan, D.W. 1930 The Theory of the Potential, Dover Publications, 
Inc. New York, 1958 

[17] Maor, L, Lahav, O. 2005, JCAP 7, 3 (ML05) 

[18] Mota, D.F., van de Bruck, C. 2004, A&A 421, 71 

[19] Peebles, P.J.E. Ratra, B. 1988, ApJ 325, 17 

[20] Peebles, P.J.E. Ratra, B. 2003, Rev. Mod. Phys. 75, 559 

[21] Percival, W.J. 2005, A&A 443, 819 

[22] Ratra, B. 1988, Peebles, P.J.E. 1988, PhRvD 37, 3406 

[23] Rubiho-Martin, J.A., Rebolo, R., Carreira, P., et al. 2003, MNRAS 341, 
1084 

33 



[24 
[25 
[26; 
[27; 

[28; 

[29 
[30 



Sievers, J.L., Bond, J.R., Cartwright, J.K., et al. 2003, ApJ 591, 599 
Spergel, D.N., Verde, L., Peiris, H.V., et al. 2003, ApJS 148, 175 
Wang, L., Steinhardt, P.J. 1988, ApJ 508, 483 
Wang, L. 2006, ApJ 640, 18 

Weinberg, N.N., Kamionkonsky M. 2003, MNRAS 341, 251 
Wetterich, C. 1988, Nucl. Phys. B 302, 668 
Wetterich, C. 1995, A&A 301, 321 



Appendix 

A. On the virial theorem for subsystems 

Let i and j denote concentric, spherical-symmetric, homogeneous subsys- 
tems, Ri and Rj, Ri < Rj, the related radii, and let the potential and the 
potential energy terms be expressed by Eqs. (2), (8), (9), (10), and (14). The 
derivation of the potential energy terms with respect to the inner or outer 
radius yields: 

dVti dVti 
Rl dR~ = ~ Qi ; Rj dR- = ° ; (82) 

Riot*: = -n,- ; ^ = o ; (83) 



Ri <)\Y,, li, OWji 1 Vij 



1 + 3wj dRi 1 + 3wi dRi 2 1 + 3w-j 



(84) 



R 3 dW 3l _ R 3 d\Y„ V %3 (h_ y2 3 ^ 



1 + 3wi dR 3 1 + 3wj dRj 1 + 3^^4 4 



dVa dVa , , 

R i*£ = 2V v ; Ri*£ = -W ij ; (86) 



dRi J ' dRj 



R >dR 3 ~ ~ X% 



5 o 15\ /5 2 3\ 1 + 3wi 
4 V -^) + {A y -a)TT^ 



(87) 



1 dRi 3 \ 2 2 1 + 3w J ' K ' 
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which, in turn, produce: 

dtt u dtt u . . . . 

dW dR~ = ' u = l >3 ' U = J^ ; ( 89 ) 

dV uv dV uv 

Ru ~dlt + Rv int = ~ Vuv ' u = i >J ; ; V = J^ ; (9°) 

owing to Eqs. (19), (89), and (90), the virial theorem for subsystems may be 
formulated as: 

_ „ dVl u „ <9fL „ dVL, „ dV uv 

r\rp p <* p « p p av n 

(9i? u <9-R„ dR u dR v 
u = i,j ; v = j,i ; (91) 

and the sum of the two alternative expressions reads: 

2r - fl -a«-- fl 'al- = ; (92a) 

T = Ti + Tj ; (92b) 
Q = Q. + + V^-j + fij- = fii + Wij + + flj ; (92c) 

where T and f2, according to Eq. (22), are the kinetic and the potential 
energy, respectively, of the whole system. 

The explicit expression of the last two terms on the left-hand side of 
Eq. (92a), and the related sum, may be calculated using Eqs. (82)-(88). The 
result is: 

r, dtt „ 1 / 1 + 3wi\ Tr 

^=-"'-(4" -iH 1+ TT3^r j ; (93b) 

„ dtt n dtt „ „ /5 2 1\ / l + 3^\ , , 

in the special case where the two subsystems fill the same volume, Ri = Rj = 
R or y = 1, Eq. (18) holds and Eqs. (92) reduce to: 
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£5, = " \ v " ~ \ v " • (94b) 

& ( fl -H) + ^ = - a - - v » - v >- - a - = - a ;(94c » 

in conclusion, the total potential energy of the whole system must be con- 
ceived as dependent on four independent variables, Vt = Q(Mi, Mj, Ri, Rj). 

On the other hand, one-component systems are subjected to no tidal 
potential, then the potential energy coincides with the self potential energy, 
which depends on two independent variables, Q = Q(M,R). Accordingly, 
Eq. (92a) translates into: 

2T-«f^0; (95) 

which cannot be splitted as a sum of different contributions, as done in 
e.g., ML05, unless it is conceived as a function of four independent variables 
instead of two, and Eq. (92a) is used instead of Eq. (95). 

B. The function f{x) 
Let us define the function: 

f( x ) = 1+x l + (l + 3w)s 

T{X) 1 + m l + (l + 3w)m ' 1 aJ 

< x < m ; < f(x) < 2 ; -1 < w < -- ; (96b) 
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where the values at the extrema of the domain are: 



1 + (2 + 6w)m + (1 + 3w)m 2 
and keeping in mind that the equation: 

{l + x)[l + (l + 3w)x) = ; (98) 
admits the real solutions: 

xi = -1 ; x 2 = — ; (99) 

l + 3w ' K J 
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the sign of the function, f(x), is determined by the inequalities: 



f(x) > 
f(x) < 
f(x) > 



< m < 



1 



1 + 3w 



< x < 



-1 



l + 3w 



m > 



-1 



m > 



1 + 3w 
-1 



< x < 

-1 

1 + 3w 



-1 



1 + 3w 

< x < m 



(100a) 
(100b) 
(100c) 



l + 3w ' 

and the function is null at x — — 1 and a; = — 1/(1 + 3w), respectively. In 
the special cases, m — — 1 and m = — 1/(1 + 3w), the function diverges 
everywhere within the domain, except at x = m where f(m) = 1. 
The first and the second derivatives are: 

df (2 + 3w) + 2(1 + 3w)x 

dx~ 1 + (2 + 3w)m + (1 + 3w)m 2 
d 2 f 2(1 + 3w) 



da; 2 1 + (2 + 3w)m + (1 + 3w)m 2 
where the first derivative is null at the extremum point: 



(101) 
(102) 



12 + 3w 



x' = 



21 + 3w ' 

and the sign of the second derivative is defined by the inequalities: 



da; 2 
d 2 / 



< 



> 



< m < 



-1 



1 + 3w 



m > 



-1 



(103) 

(104a) 
(104b) 



da; 2 1 + 3w 

in the range of interest, — 1 < w < —1/3. Accordingly, the extremum point, 
x\ is a maximum and a minimum, respectively. 

At the extremum point, the function attains the value: 

Q in 2 1 
/(* f ) = -T i o ; ■ ,o ■ o ^ ■ ,n ■ o ^ 5 5 (105) 



4 1 + 3w 1 + (2 + 3w)m + (1 + 3w)m 2 

on the other hand, the extremum point is attained at the upper limit of the 
domain, x^ = m, provided the parameter, w, has the value: 

2 m + 1 



w = — - 



3 2m+ 1 



(106) 
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in particular, w = —2/3 for x* = m = and w — —1/3 for x* = to — > +oo. 
Within the range, — 1 < w < —2/3, the extremum point lies outside the 
domain, x^ < 0. 

To gain further insight, a few special cases shall be studied with more 
detail, namely w = —1, —2/3, —1/2, —1/3, the last to be conceived as an 
interesting limiting situation. The general case, —1 < w < —1/3, is expected 
to show similar properties as in the closest among the above mentioned spe- 
cial cases. 

In the special case, w = —1/3, Eq. (96a) reduces to: 

f( x ) = ; (107) 

JK ' 1 + m K J 

which is positive within the domain, < x < m, according to Eq. (100a). The 
extremum point is a maximum and occurs at infinite, according to Eqs. (104a) 
and (105). Then it belongs to the domain only in the limit, to — > +oo. The 
function, f(x), for different values of the parameter, to, is represented in 
Fig. 7, top left. 

In the special case, w = —1/2, Eq. (96a) reduces to: 

f(x) = 9 2 ± X ~ X \ ; (108) 

2 + TO — TO 

which, according to Eqs. (100), is positive within the domain, < x < to, for 
to < 2, and within the domain, 2 < x < to, for to > 2; on the other hand, 
it is negative within the domain, < x < 2, for to > 2. The (allowed) zero 
of the function occurs at x — 2. In the special case, to = 2, the function 
diverges everywhere within the domain, except at x — 2, where /(2) = 1. 

Owing to Eq. (103), the extremum point occurs at x* — 1/2, which is a 
maximum for < to < 2 and a minimum for to > 2, according to Eqs. (104). 
The value of the function at the extremum point is: 

/W = ~ — 2 ^ * t = ^ ; ( 109 ) 

4 2 + to — m 2 2 

which belongs to the domain only if to > 1/2, to 7^ 2. The function, /(x), 
for different values of the parameter, to, is represented in Fig. 7, top right. 
In the special case, w = —2/3, Eq. (96a) reduces to: 

tt*) = \ z ^* ; (no) 

1 — rrr 
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Figure 7: The function, f(x), for different values of the parameter, m, within 
the domain, < x < m, with regard to special choices of the parameter, w, 
(from top left in the clockwise sense): —1/3, —1/2, —2/3, —1. The value of 
m related to each curve, starting from the left, is m = i/4, < % < 12. The 
horisontal axis corresponds to m — > +oo. The locus of ending points is the 
horisontal line, f(x) = 1. In the special case, m — 0, the curve reduces to 
a single point, (0, 1). The locus of extremum points is marked by a dashed 
vertical line for w = —1/2, and coincides with the vertical axis for w = —2/3. 
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which, according to Eq. (100), is positive within the domain, < x < to, for 
to < 1, and within the domain, 1 < x < to, for to > 1; on the other hand, 
it is negative within the domain, < x < 1, for m > 1. The (allowed) zero 
of the function occurs at x = 1. In the special case, to — 1, the function 
diverges everywhere within the domain, except at x — 1, where /(l) = 1. 

Owing to Eq. (103), the extremum point occurs at x' = 0, which is a 
maximum for < m < 1 and a minimum for m > 1, according to Eqs. (104). 
The value of the function at the extremum point is: 

/(xt) = — ; xt = ; (111) 

which, in any case, belongs to the domain. The function, f(x), for different 
values of the parameter, to, is represented in Fig. 7, bottom right. 
In the special case, w = —1, Eq. (96a) reduces to: 

/(*) = / = x - l x \ ; (H2) 

1 — to — 2m 2 

which, according to Eqs. (100), is positive within the domain, < x < m, 
for m < 1/2, and within the domain, 1/2 < x < m, for m > 1/2; on the 
other hand, it is negative within the domain, < x < 1/2, for m > 1/2. 
The (allowed) zero of the function occurs at x — 1/2. In the special case, 
m = 1/2, the function diverges everywhere within the domain, except at 
x = 1/2, where /(1/2) = I. 

Owing to Eq. (103), the extremum point occurs at x^ = —1/4, which is 
outside the domain. The function, f(x), for different values of the parameter, 
to, is represented in Fig. 7, bottom left. 

C. The effect of partially clustered quintessence on the 
virialized configuration 

With regard to a virialized configuration in the special case of fully clustered 
quintessence, let be the energy change, E — E FC , which defines the viri- 
alized configuration in the case of partially clustered quintessence, related to 
A/j, = CVfc"^ according to Eqs. (54)-(57), (65), and (69). The combination 
of Eqs. (56), (65), and (69) yields: 



An = - 



3GMI 



rn 



5 Rmax Vfc 



[1 + 3w)C(2 - Qx 1 + (2 + 3w)(x ; (113a) 
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x = mr] F Q W ; < x < m ; (113b) 
and the substitution of Eq. (113b) into (70) produces: 

<f>(x; w, C) = (1 + Sw)((2 - Qx 2 + (2 + 3w)(x ; (114) 
which makes Eq. (113a) be cast into the form: 

Afi 3^, ; (115) 

accordingly, the sign of the energy change, AQ, is opposite to the sign of 
the function, (j>(x). Keeping in mind that positive and negative AQ imply 
expansion and contraction, respectively, with regard to the virialized config- 
uration in the special case of fully clustered quintessence, the effect of partial 
clustering may be deduced from the sign of <f>(x). According to Eq. (114), the 
solutions of 4>(x) = are: 

1 2 + 3m , , 

xi — ; x 2 = ; (116) 

2-Cl + 3w ' v ; 

and the sign of <f>(x) is defined as: 

(f>(x) > ; mm(xi,x 2 ) < x < max(xi, x 2 ) ; (117a) 
(f)(x) < ; x < min(a:i, x 2 ) ; x > max(xi, x 2 ) ; (H7b) 

where < x < m in the case under discussion, conform to Eq. (96). 
The first derivative: 

-^ = 2C(2-C)(l + 3w):r + C(2 + 3 W ) ; (118) 

is null at the abscissa: 

X 22-Cl + 3^ ' [ yj 

where, owing to Eq. (116), — Xi\ = \x^ — x 2 \, or x^ = x 2 /2. 
The second derivative: 

-0 = 2C(2-C)(l + 3«;) ; (120) 
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is negative within the range of interest, — 1 < w < —1/3, < ( < 1, which 
implies that the extremum point is a maximum. The substitution of Eq. (119) 
into (114) allows the calculation of the function at the maximum. The result 
is: 

, (l t ) = _i_CJW . (121) 

vy ! 42-C l + 3w ' y ! 

which is non negative, 4>{x^) > 0, in the case under discussion. 
Further inspection of Eq. (119) shows that: 

x*>0 ; -\ <w< -\ i ( 122a ) 
< ; -1< w < -- ; (122b) 

x* = ; w = ~\ 5 ( 122c ) 
xt - +oo ; w - ; (122d) 

independent of the value of £ within the assigned range. 

The particularization of Eqs. (114), (116), (119), and (121), to a few spe- 
cial values of (, yields: 





4>{x) -- 


= (x 


1 




0(xt) 


— > +oo ; x' = 


in the special 


case, w = 


-1/3; 




<f>{x) = - 




(2 - C)x 2 - x 




0(xt) = 


1 c t 

82-C ' X ~ 


in the special 


case, w = 


-1/2; 



(123a) 

+oo ; (123b) 



; (124a) 



0(x) = -C(2-C)^ 2 ; (125a) 
0(x f ) = ; ^ = ^2 = 0; (125b) 
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in the special case, w = 



-2/3; 



4>{x) = -C [2(2 - ()x 2 + 



x 



(126a) 




42-C 



1 1 



(126b) 



in the special case, w — — 1. 

The function, 0(x), expressed by Eqs. (123a)-(126a) and related to dif- 
ferent values of (, are represented in Fig. 8 with regard to the special cases 
considered above. 

Keeping in mind Eq. (115), the transition from a virialized configura- 
tion related to fully clustered quintessence to its counterpart related to par- 
tially clustered quintessence, implies expansion (0 < 0) for values of the 
quintessence equation of state parameter, w, in the range, — 1 < w < — 2/3. 
On the other hand, both expansion and contraction (0 > 0) may occur for 
—2/3 < w < —1/3, according if the variable, x, is sufficiently distant from 
and/or the parameter, w, is sufficiently close to —2/3, and vice versa. 
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Figure 8: The function, <f>(x), for different values of the parameter, (, with 
regard to special choices of the parameter, w, (from top left in the clock- 
wise sense): —1/3, —1/2, —2/3, —1. Dotted curves correspond to £ = 1/4, 
dashed curves to ( — 1/2, dot-dashed curves to ( — 3/4, and long-dashed 
curves to ( — 1, which is the limit of fully de-clustered quintessence. In 
the limit of fully clustered quintessence, £ = 0, all curves coincide with the 
horisontal axis. 
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